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APPARATUS AND METHOD FOR DETECTION OF ONE LUNG INTUB A l ION 

BY MONITORING LUNG SOUNDS 

FIELD OF THE INVENTION 

The present invention relates to acoustic detection of one lung intubation in ventilated patients. 

BACKGROUND OF THE INVENTION 

During general anesthesia, for proper air way management, an endotracheal tube is inserted 
into the patient's trachea through which the patient is ventilated. The tube is inserted during the 
primary induction and placed so that its tip is located above the carina - the bifurcation of trachea 
into the two main bronchi. The location of the tip of tube is critical: it should be placed, and 
maintained above the bifurcation. A correct position of the tube, in which both lungs are ventUated, 
is called Tracheal Intubation (TRI). If the tube is misplaced or shifted due to patient movements, 
cases of One Lung Intubation (OLI) may occur. Prolonged cases of OLI should be avoided since .t 
may cause insufficient oxygenation and may damage the non-ventilated lung. OLI was found to be a 
cause of desaturation and a cause of malfunction during anesthesia, and there is currently no reliable 

device or method for detecting OLI situations. 

Currently known methods for detecting one lung intubation including the stethoscope and 
capnograph, have proven either unreliable. Pulse oximetry is the most reliable known method but 
provides results with latency of 2 to 5 minutes, which may be too long to prevent damage. There » 
an ongoing medical need for methods and devices for detecting one lung intubation in real time. 

The following published documents provide potentially relevant background art and are 

incorporated herein by reference: 

Sod-Moriah G„ Gelber O, Gurman G. and Cohen A. "Monitoring of Separate Lung 
Ventilation in Anesthesia and Intensive Care" Proc. of the IEEE 18th convention of 
electrical and electronics engineers in Israel, Tel-Aviv, March 7-8, 1995. 
Sod-Moriah G., Cohen A. and Gurman G„ "Detection of One Lung Intubation 
Incidents in General Anesthesia and Intensive Care," Proc. of the 13th Int Conf. 
BIOSIGNAL 96, pp. 282-284, Brno, Czech Republic, 1996; 



WO 2005/094179 



PCT/EL2005/000369 



2 

US 2003/0018276 of Mansy H. and Sandler R, titled "Acoustic detection of endotracheal 
tube location." 

It is also noted, that a poster presentation Gurrnan et al t "Continuous Monitoring of Separate 
Lung Ventilation" ASA annual Meeting Orlando Florida October 12-16 2002 provides potentially 

5 relevant background art- 
US 2003/0018276 of Mansy HL and Sandler R. titled "Acoustic detection of endotracheal 
tube location" purportedly disclosed a system 2nd method for use in detecting an endotracheal rube 
location within a body by electronically detecting breath sounds. Unfortunately, insufficient details 
were provided to allow others to reproduce the disclosed results. 

10 There is an ongoing medical need for methods and apparatus for reliably detecting a one lung 

intubation condition in ventilated patients. Preferably, the method and apparatus will be operative in 
the presence of background noise of an operating room or intensive care ward. 

SUMMARY OF THE INVENTION 

] 5 it is now disclosed for the first time a method of detecting a one lung ventilation situation in 

a human subject. This method includes detecting indigenous lung sounds emanating from a region of 
the body with acoustic sensors to produce an electronic signal, and generating an output indicative of 
the one lung ventilation situation by processing said detected indigenous lung sounds. 

According to some embodiments, the processing includes computing an autoregressive 

20 moving average (ARMA) or autoregressive model of the electronic signal. 

Not wishing to be bound by theory, it is understood that the human body is not composed of 
a uniform medium, but is heterogeneous. Local acoustic properties vary between different types of 
tissue. As such, noise generated by sources such as the lungs is subjected to a certain amount of 
dispersion as the noise is transmitted through the human body. Parts of a specific noise generated by 

25 a source and transmitted through the body thus reach a detector on the surface of the human body at 
different times. Therefore, it is now disclosed that there is a correlation between a measured noise 
signal emanating from the lungs and the history of the measured noise signal. 

As such, non-linear models or even linear models such as autoregressive moving average 
(ARMA) models or autoregressive models provide a reasonable representation of a source noise 

30 signal, and are useful for determining whether detecting lung sounds are from one or two intubated 
lungs. 
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In some embodiments, the disclosed method includes determining a number of active 
distributed noise sources of indigenous lung sounds or a number of distributed random sources m 
order to detect a one lung intubation situation, wherein a detection of one active distributed no.se 
source is indicative of OLI while two active distributed noise sources indicates TRL In some 
embodiments, the detection is carried out by the general approach of Blind Source Separation 
Unlike previously disclosed methods for Blind Source Separation, the currently disclosed method 
does not require a determining of the actual signal generated by the distributed noise source, 
Furthermore, the methods of the present invention are appropriate for any distributed noise source. 

Furthermore, it is noted that the present invention is not limited by the aforement.oned 
techniques, and that any method known in the art for obtaining or recovering any or part of spat.al 
statistics of distributed noise sources from measured lung noises is appropriate for the present 
invention. In some embodiments, the processing of lung sounds includes using neural networks m 
order to obtain spatial statistics of the of the lung noises. In some embodiments, the processing of 
iun* sounds includes linear or non-linear modeling of the lung sounds such as Hidden Markov 
Model (HMM). In some embodiments, the processing of lung sounds includes blind deconvolufon 
and system identification using higher-order statistics. 

It is now disclosed for the first time a method of detecting a one lung ventilation .human m 
a human subject. The presently disclosed method includes electronically detecting indigenous lung 
sounds emanating from a region of the body, and generating an output indicative of the one lung 
ventilation situation by processing the detected lung sounds. 

According to some embodiments, the detecting includes receiving a plurality of electncal 
signals from a plurality of acoustic sensors, and at least one acoustic sensor is disposed adjacent to at 
least one region selected from the group consisting of a chest region of the body and a back regjon of 
the body. Exemplary locations for the acoustic sensors include the left side of the chest, the nght 
side of the chest, the left side of the back and the right side of the back. 

In some embodiments, the processing includes processing only electrical signals from 
acoustic sensors placed on the back. In some embodiments, the processing includes processing only 
electrical signals from acoustic sensors placed on the chest. 

As noted earlier, local acoustic properties vary between different types of tissue leading to 
multipath propagation of lung sounds. The detected sounds can be described as a convolve 
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mixture and modeled by a convolutive mode!. As such, the detected signal can be said to have 
memory. 

In some embodiments, the detecting includes receiving a plurality of electrical signals from a 
plurality of acoustic sensors, and the processing includes computing a parameter indicative of a 
5 relation between a received electrical signal and a past and/or future behavior of the rer^A 
electrical signal. For embodiments where a different electrical signal is received from each acoustic 
sensor, it is noted this parameter may itlatc any of received electric signals with the past and/or 
future behavior of the same or any other received electric signal- 
In some embodiments, the processing includes computing a parameter indicative of a relation 
10 between a received sienal and at least one of a history and * future behavior of the received 5™!. 

In some embodiments, the parameter is indicative of a relation between an electrical signal 
from a Grst acoustic sensor among the plurality of acoustic sensors, and the history and/or future 
behavior of the electric signal from the same first acoustic sensor. 

In some embodiments, the processing includes computing a parameter indicative of a relation 
15 between a received electrical signal during a first time window and a received electrical signal 
during a second time window, wherein the first and second time windows may overlap. 

In some embodiments, the first and second time windows overlap by at least 0.5 seconds. 
Alternately, the first and second time windows overlap by at least 0.75 second. Alternately, the first 
and second time windows overlap by at least 1 second. Alternately, the first and second time 
20 windows overlap by at least 1.5 seconds. Alternately, the first and second time windows overlap by 
at least 2 seconds. 

It is noted that the time window overlap is an optional feature, and in some embodiments 
there is no overlap whatsoever between analyzed time windows. 

In some embodiments, the parameter indicative of a relation between a received electrical 
25 signal and a history of a received electrical signal is computed for a plurality of times. 

In some embodiments, the relation is indicative of a conditional probability relation such as a 
conditional probability that a future signal will have a certain form given the form of the present 
signal. 

In some embodiments, the processing includes computing a parameter related to a covariance 
30 matrix of said conditional probability relation. In one particular embodiment, the covariance matrix 
is a residual covariance matrix. 
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One exemplary parameter indicative of a one lung ventilation situation is the magnitude of an 
M, of the residua! covarianee matrix. Selection of the specific eigenvalue for detecting one 
tag intubation depends on the speoific number of acoustic sensors employed. In some 
embodiments, a lower value of*, eigenvalue of the residual covarianee mart* is mdtcatiue of the 



one lung ventilation situation. 



obtained. 



g ventilation situation. 

Although certain embodiments of the present invention relate to techniques of Blind Source 
Separation, recovery of original transmitted lung sounds is not a requirement of the present 
invention. According to some embodiments, the processing includes determining only a number of 
distributed noise sources of said indigenous lung sounds. 

Nevertheless, in some embodiments an estimate of the originally transmitted lung sounds « 

According to some embodiments, the processing includes obtaining a parameter indicative of 

spatial statistics of the indigenous lung sounds. 

In some embodiments, the spatial and/or temporal statistics of the indigenous lung sounds ,s 
monitored in time, and deviations in the spatial statistics of the lung sound are indicative of a change 
Ln intubation status, e.g. a change from TRI to OIL In one particular example, the spatial stat.st.cs of 
the lung sounds are obtained during a time period of known TRI such as at the beginning of surgery 
to train the system. At a later time for which it is desired to detect the presence or absence of the one 
lBagJMBti | dfiQ Situation Statistics ofthe Inn, sounds are com P ared_wit^^ 

the, time per ^H of known TRI. 

According to some embodiments, the processing includes determining a source scattenng 
parameter indicative of a scattering (spatial distribution) of noise sources. More scattered, non- 
coherent point noise sources are more indicative of TRI, while less scattered noise sources are 
evident by a smaller second eigenvalue of the residual covariance matrix. 

In experiments conducted by the present inventors, a one lung ventilation situat.cn m a 
human subject was determined in the presence of uncancelled, random background noise associated 
with an operating room or intensive care ward. Although certain types of background noise can be 
cancelled out using techniques such as simple filtering or adaptive noise cancellanon, and the 
presently disclosed methods and devices do not preclude usage of these techniques, not every 
random background noise is necessarily cancelled out, especially if this noise is not treated using the 
aforementioned techniques. The present invention provides methods and devices for detecting a one 
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lung ventilation situation even in the presence of uncancelled, random background noise of a 
loudness associated with an operating room or intensive care ward. 

In some embodiments, the processing includes processing the detected indigenous lung 
sounds in a way that is insensitive to uncancelled, random background noise of a loudness associated 

with an operating room. 

In some embodiments, uncancelled random background noise includes at least 70 decibels of 
uncancel noise, or at least 75 decibels of uncancelled noise, or at least SO decibels of uncancelled 



noise. 



Optionally, the stage of detecting includes detecting noise other than lung sounds, and the 
stape of proving includes using an adaptive filtering technique to filter noise. 

In some embodiments, a specific parameter indicative of a one lung intubation situation is 
calculated, and when the value of the calculated parameter drops below or climbs above the 
threshold of a predetermined threshold value, output indicative of a one lung intubation situation is 
generated. Thus, by adjusting this predetermined threshold value, it is possible to configure the 
device such that an occurrence rate of false alarms or false positives and missed OLI situations or 
false negatives changes according an error tradeoff curve. 

For example, in some embodiments, a smaller value of the calculated parameter is indicative 
of a one lung intubation situation. Thus, by raising the predetermined threshold value, instances of 
one lung intubation are less likely to be missed, and the occurrence rate of false negatives drops, 
while concomitantly, the device is more likely to generate false alarms or false positives. 

Thus, in some embodiments, the processing unit is adapted such that at most 9% of 
identifications of OLI are false positive identifications, and at most 2% of said identifications are 

false negative identifications. 

Alternately, the processing unit is adapted such that at most 4.5% of identifications of OLI 
are false positive identifications, and at most 4.5% of said identifications are false negative 
identifications. 

It is now disclosed for the first time a method including selecting a population of human 
subjects sufficiently large to give statistically significant results, and identifying a one lung 
intubation situation in a subpopulation of the population, wherein at most 9.6% of the identifications 
are misidentifications. 
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„ is now disclosed for ft. tat time a medtod including selecing a population of human 
. " I"! lar „ e „ give statistically significant results, and IdeoUfving a one lung 

:i vr—ation ^ ^ - -. - - r - -~ 

„ false positive identifications, and a, most 4.8% of the tdenettlcattons are fa.se ne a at, 

""Til now disclosed for the firs, time a meftod including selecting a popu.ation of huntan 

subj J fficiently large to give statistic* si^ean, results, 

J bali o„ si.ua.ion in a subpopuladon of the popuiation, wherein a. nros, M of fte 

„ fa,se positive iden.Mea.ions, and a. nros. 2% of ft. iden«~ are false neganve 

embodiments, a popuiadon sufucienu, large .o give — 
teiudes a. least 20 individuals. A,.en,ate,y, the popuiadon includes at leas 50 individual^ 
mJI,,. d,e populadon includes - leas, 200 individual Alternately, .he popu.ar.on .ncludes a. 

~ ' nl d "tl embodiments wii, he appare, fro. fte detailed descripdon and espies 



that follow. 




FIG. 1 provides block diagram of an exemplary MIMO AR model. 

FIG. 2 provides a graph of P t of GLRT for coherent sources. 

FIG 3 orovides a graph of eigenvalues of R versus the scattering level, «. 

HO 4 ^ovides a sLnatfe of exempt locutions of microphones on the hack of me patten. 
Lordin. to some embodiment of me present invention. The TRi si.ua.ion is iliustrated m F G. . 
^. 5 ""us.ra.er some record* breathing cycles bods for bod, one lung inbrba.ion and trachea, 

intubation cases. ^ 

FIG . 6 provides an exemplary graph of the second largest eigenvalue of R fcnct.cn of .me 

both one lung intubation and tracheal intubation cases. 

FIG.7 provides the DET of a classifier based on the second highest eigenvalue of estimated R . 
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™ ggge events of the present inven.cn that 
computing certain autoregression functions of a detected acoustic signal enables detection of a one 
lung intubation situation, even in the presence of background noise associated with operating rooms 
and intensive care wards. In particular, an algorithm for detecting the number of ventilated lungs 
from recorded breathing sounds has been developed. In some embodiments, this algorithmjtfsumes 
* MIMO (Multiple Input Multiple Output) system, in which a m»!t!-d?mens5r,r, a s AR (Auto- 
Regressive) model relates the input (lungs) and the output (recorded sounds). The unknown AR 
parameters are estimated, and a detector based on the estimated eigenvalues of the residual 
covariance matrix is developed, in order to detect a one lung ventilation situation. 

In the examples presented, a number of noise sources is estimated using measures such as the 
Akaike Information Criterion (AIC) or the Minimum Description Length (KiDL). All of these 
measures make several assumptions about the sources. Unfortunately, the problem at hand does not 
obey these assumptions, the most notorious of which is the assumption of coherent distributed noise 
source. The lung is a diffused source rather than a point source. 

Under the assumption of coherent distributed sources the large eigenvalues of the residual 
covariance matrix correspond to the sources and the small ones correspond to the diffuse noise. In 
the case of diffuse (distributed), non-coherent sources, the distinction between these two groups is 
not clear. Threshold methods are used in order to estimate the number of sources. 

Although in some embodiments this algorithm derives from a Blind Source Separation 
model, a presently-disclosed algorithm estimates only the number of active sources or lungs, and 
does not require estimation of the source signal itself. The source signal is transmitted via the chest 
or back of the patient to the sensor. Although there are non-linearities associated with transmission 
channel, it has been discovered by the present inventors that assuming a linear transmission channel 
is functional for detecting one lung intubation. Furthermore, although the most general discretion of 
a linear channel is an ARMA (Auto Regressive Moving Average) with poles and zeros, it is known 
that a less general high-order AR (Auto Regressive) model including only poles provides an 
approximation to an ARMA model. In accordance with some embodiments of the present invention, 
it is disclosed that the AR model is functional for detecting one lung intubation. 

The methods of the present invention are appropriate for any distributed noise source. In 
some embodiments, a "distributed noise source" as used herein is composed of point noise sources 
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f t la* 75 cm 2 In some embodiments, a "distributed noise source" as used 
distributed in an area of at least 7. cm - in som 

. , • t nnit . e sources distributed m an area of at least _uu 

Exiles of computed spatia, include but a« no, limited to a cross/jom. covananco 

of the measurement processes at the different sensors. „„,;„«;„„!„ 

The followin* examples are to be considered merely as M» and non-lrmmng 
nalure r^LTien, to one sHlled in the en to which the present invenuon pertains tha, many 
Z£L. peruLiona, and variations may he made wirhou, departing from the scope oftbe 



invention. 

EXAMPLES 



pf^fPfF' MOOF.T, FORMULATION 
h me present example, the breathing sound signals are recorded by 4 microphones artached «o 
0,0 £j, hac, Previous atempts » derec, 0U by comparing rhe amphhrde o *e reco^ 
sou 1 in right and left sides did no, resu,, in reliable methods, because each one of *e mrcroph nes 

sounds generated by horn ,un g , .„ order ,o overcome mis P ro„em, a convolunve - 
"broach is present. In me cunan, examples, an AR mode, met relates rhe lungs and Ore 

"™ - - — « *- * - - m * use : : p ' ica,i :i 

^ and audio processing and its computation, complexity is relatively simp*. k d» mode, 
Ite«i,a,ed lu , represents a source. Our goal is to detec, a situation of which only one lung ,s 
Elated, from the received signals by ,he sensors. « is assumed the, rhe 
ventilated lungs ara independent. Fig. , shows a blocx diagratn of *e proposed MD£»-»* • 
in which m represents the sources (lungs), and yW represent me sensor (mrccophones) 



measurements. 
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Let K and L denote the number of sources (lungs) and sensors (microphones), respectively 
(K<L). Therefore, the vector of source signals, x[w], is defined as a K x 1 vector as follows: 

= [*,[»] -^["l • - • *<[»]] . 



(1) 



Thelxl measurement vector is defined as: 



yM -[*[»] y 2 M ■ • ■ nWf 

The relation between the source signals and the measurements is assumed to be given by a MIMO- 
AR model: 



y[/7] = Ay (w) [n] + CxM + e[/j], 



where y 'M is an M,x 1 vector defined as follows: 



(4) 



and y# W W is an Wxl vector which contains the past values of the i-th sensor, y,{ri\ 9 up to sample 



M: 



(5) 



A is an L*A4L matrix defined as: 



a 



a 



(6) 



where a# is an Mxl vector, which relates the samples of the Mh sensor, with the past values 
of the y-th sensor, , C is an Lx K matrix whose y-th element relates the samples 

of source j and sensor /. Finally, e[n] is an Lxl vector representing additive white noise. It is 
assumed that the noise and source signals are independent, zero-mean, Gaussian with covariance 

matrices <x I and I, respectively. The last assumption can be employed with no loss of generality, 
because the covariance of the sources is determined by the matrix C, as it can clearly be seen from 
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, ♦ -u f-«« yMIv^M is Gaussian: 
(5) . As a result, it is obtained chat the conditional dilution of 

y[«31y (W) [»] ~ ^(Ay (W) [»l» R )' where R is defined as: 

R = CC r +cr J I. < 7) 

^meters- A, R U and K must be estimated from a set of 
It is noted that the unknown parameters. A, k, m 

r,l vfM It is also assumed that all the initial conditions are zero, ,e. *W»H 
5 measurements, y [1], • -,y[W- » « 35 SUCC essful estimation of K, the 

for n<0, and that the input and noise signals are stannary . In fact, successful 

number of sources (lungs), is the key for the OLI detection. 

•gVATVTPT.E 2. TTTP M0[ ESTIMATOR 

r ~c r we need first to estimate the unknown matrices, 
b order to determine the number of sou *s, K we need fa Max5mum _ Likelihood 

10 A and R, from the //samples of the data: jPl^-jW F ^ > ^ 

nmr n erimatnr is used The ML estimator of the matnces A and R, is obtamea oy 

(ML) est.mator is used, me ^ giyen me 

boarithm of the conditional probability denstty function u> ; 
unknown matrices, which is : 
,o g /MH y[W )|R,A) = -flog(2T)-fl°8W- 

12 .if f( yM -/^->w) r R- , (ri»i-*»*"w)] • (8) 

• hv onuatine its derivatives with respect to A and R, 
The loa-likelihood function can be maximized by equating its oer 

The log mceu Seg A dix A): 

and solving the two resulting matnx equations. This process ytei 



20 



£yWu r [»]Yt u W« r [»]] (9a) 



and 



R M =*ty[* r w- A «(*|" Wyrw ) ( 9b) 



25 



™. u,of .ode, orfe, se.ac.ion m ed,ods based on iofonnadon M cn.ena „, H .«) «-» 
r e — -W in orde, , es t in,a,e *e .node, orfe,, M and *. ~« * 
TO, me *od was deve,oped and X during our wo*. U did „.« show a «t «-» Wen 
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applied to real breathing sound signals . Therefore, a Generalized Likelihood Ratio Test based 
method was developed and tested as shown in the next Section. 



10 



EXAMPLE 3: GENERALIZED LIKELIHOOD RATIO TEST 

In the private case of lungs as sources, the number of sources can be only one or two. Therefore, 
for the purpose of decision of between TRJ case and OLI case, the GLUT is used [15]. This test 
is based on the ratio between the probability density function under each hypothesis, while the 
maximum likelihood estimator is used to estimate the unknown parameters under each 
hvpothesis. Let us denote the following? hypothesis: 

^ : Only one source exists for the system (OLI case, K=\) 

&2 : There are two sources for the system (TRJ case, K=2) 



15 



20 



The development of the Log-likelihood function under the i-th hypothesis, leads to the following 
expression (assuming the noise variance, cr 2 1 is known) : 

NL 



ML N ( rr / 5 \ l ' K \ 

log/Cy[l],.- > y[^]|R»A;//,) = -~log(2^)-ylog^n / '( cr ) J- 



(10) 



where Mm are the two highest eigenvalues of R ( 7 i and L is the number of sensors (Proof: 
See Appendix B). 

As a result, the GLRT for decision between H \ and H i is as follows: 



log 



/(yPL- ..>y[Ar]lR,A;/y,) 



.1-2 \ 



/(m yt«]IR.A;# s )J 2 [ A("') J 



do 



As can clearly be seen from (11), the second highest eigenvalues of R is actually the detector of 
OLI situation, under the assumption that or 1 is known, and that the sources are point sources. 
Simulation results given in the next section show the performance of the proposed detector under 
coherent and incoherent distributed sources assumption, while the last assumption is a more accurate 
25 model for lungs sources. 
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EXAMPLE 4: SIMULATION RESULTS 
/. Coherent sources f . . 

h order ,o evaiuate the perfbnnence of *e estimators as a taction of the p— of *. 
TO del simu>ations with synthenc dam were ^ 71- MIMO-AR syaem as MM . (3). 

1,- and the ML — - of R and A were Ccuiated - <»• - * — 

*e p _ of the system were as foHows (unless otherwise Is indmated): the number of sources 
was K-l, me number of sensors was i-4 with AR order, M-5 and noise variance, <r -l.Tbe 
matrix C was chosen to be: 

0.3 0.15 



C = 



0.15 0.3 
0.9 0.6 
0.6 0.9 



mi *. nratri* A was chosen to be stable representing system poles inside the unit circle. Simulation 
results of the ML estimators of the matrices A and R, can be found tn - 

The behavior of the GLRT as a taction of the number of independent samples, M . -mm*. 
Tha second highest eigenvalue, <>, was extracted and compared to a threshold value of the noise 
,evel, *' . A total number of .M000 iterations for each N were performed, and the probabthty of 
error of K is defined as: 

number of uncor rectly estimated K 

p <= T~ 0 2 > 

F5g 2 shows the probability of error, , as function of number of samples, N. It can be seen that 
probability of error decrease as the number of samples grows. This fact justifies the use of a 
threshold value of the noise level when the sources are coherent. 

2 Incoherent Distributed Sources Simulation Results 

It is well known that each lung is composed of several independent point sources, and 
therefore should be treated as distributed sources. In order to evaluate the performance of the 
al corithm under this assumption, spatially distributed sources were synthesized. In th,s 
simulation, two distributed sources were synthesized. Each distributed source was composed of 
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four independent point sources with close spatial signatures. Therefore, the columns of C were 
chosen to be: 



C = 



c, c,+£Ac„ c, + 2sAc, 2 c,+3gAc„ c 3 c,+26Ac,, c^fAc^ cyfreAcg 



source I 



J(13) 



where c i and c * are orthogonal, and i c < Ac * Ac " Ac ' 3 ^ is an orthonormal ?rou». e is » 
5 constant which determines the distribution width. The vector of sources x[n) represents eight 
independent sources. Therefore, the product Cs M represents two distributed sources, where each 
one is constructed by 4 independent sources. 

In order to examine the performance of the GLRT under incoherent distributed sources 

assumption, a comparison between the eigenvalues of R for cases of one and two distributed 

10 sources is given. In Fig. 3 the eigenvalues of R are drawn as function of £ for K=l and K=2. It can 
be seen from Fig. 3 that even when the sources are widely distributed, the eigenvalues of the source 
signals sub-space are separated from the eigenvalues of the noise sub-space. Therefore, despite the 
fact that the lungs do not function as a coherent sources model, it was decided to examine the second 
highest eigenvalue of R as a detector for OLI situation in real breathing sound signals. In addition, 

15 it can be also be seen the threshold value in the case where the sources are widely distributed should 

be higher than the noise level, <? =1. 

EXAMPLE 5: EXPERIMENTAL RESULTS 
In order to examine the disclosed model for OLI detection, a database of recorded breathings 
20 was established. The database was composed of 24 patients which were recorded in a surgery 
room in both situations: during correct ventilation, when the tip of the tube is placed above the 
carina, and during a situation of OLI when the tip of the tube is under the carina and only one 
lung is ventilated. 

During each experiment, the microphones were attached to the patienfs back, as shown in 
25 Fig. 4, recorded the breathing sounds of the patients in both situations. The ventilations were 
performed manually and not mechanically, in order to achieve higher signal-to-noise ratio in the 
recorded sounds, and the real position of the tube was validated each time by fiber-optic. The 
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15 



rf ,A in the main surety room of medical center Soroka - Israel, during 
experiments were performed in the main surgery 

the anesthesia part in the beginning of the surgeries. 

* order to attenuate some of the irregular background noises outside of the spectral range for 
STZ s such as monitor beeps, doctors' discussion, etc. the recorded signals were band- 

■ B— * ^ ^ b - dW5dth " 100HZ - 6 ° 0H2 - ^ ^ ^ 
ir 0 f the cut-off frequency of 600Hz, do„Ung operation with a factor of 0.3 

WM ;?sC' amplitude from each microphone depends on the particular ,ocation of the 
J£ 71 bod, the anatomy of the particular patient, and on the gain of the sampling 

^rrirefore, it turned out that each channel, output had a different signal ampli.de, and the 
irilce, I', also differed between microphone, In order to treat this problem, a 
normalization ofeach channel according to the noise level on*e channel was done 

It is also noted that the aforementioned techniques were only sufficient to reduce some 

,„h the algorithm itself was robust enough to 
ambient noise associated with an operating room, and the algorithm 

determine a OLI situation in a manner that was insensitive to irregular no.se of an operatm, 
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25 



room. 



n. «co,ded bre a*i„ g co^in bo* of OLI and TO T„e <^ 

JLed .o ou,-off » of 4kHz, and .he da. w K e divldad ioto wndows o f 2 000 ^ 
were nmiie hj hegt avai]ab!e AR 

each with 80% overlap. Because of fact that AIC and MDL nave 
:l order when applied to real data, an arbitrary AR order of 15 was set 
computation complexity and the available processing time. The unknown matnces A and R 
estimated for each window, using the ML estimator developed in Section IL 

Pig 5 shows a few breaming cycles of both OLI and TRI situations, recorded by *e fcu 
^phones after pre-processing. As it can be seen from this figure, determination between OLI and 
TRI cases by only the amplitude of the recorded sounds is not a simple task. 

, . . • . . opnvahie 0 f r as a function of time, as a result ot 
Fig. 6 shows the second highest eigenvalue oi K Mt <*> 

• r:„ < Ac it can clearly be seen from Fig. 6, OLI and TRI 
processing the measurements shown in Fig. 5. As it can cieariy 
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cases can clearly be discriminated, by the second highest eigenvalue of in every breathing 
cycle- The results of the proposed algorithm were consistent over the 24 experiments. 

Estimation of the performance of the system was performed using "Leave some out method" as 
follows. Twenty different experiments were used to extract the histograms in order to train the 
5 system. The rest of the 4 experiments were used to validate the system and were tested according to 

the extracted statistic* in iiie iiwiuiiK fnovess. This process was repeat*! q iin.es, each time a 

different' group of four validation experiment was used. As a result, a validation was performed 
using a total number of 24 experiments, in a "patient independent" mode. 

There are two types of errors in OLI detection: P miss , the probability of a true OLI to be wrongly 

1 0 detected as TRI, and P FA> the probability of TRI to be detected as OLI. The Detection Error Tradeoff 
(DET) curve is a common mean to display these errors. The DET curve provides information about 
the device's performance, where each point on the curve shows the Pfa and P miss for a given 
threshold. The threshold of a real monitoring system should be calculated according to the requested 
sensitivity of the system, while taking into consideration the allowed P miss of the system. Fig. 7 

15 shows the DET curve of the proposed decision system, which was computed according to the 6 
iterations described above. The Equal Error Rate (EER) point is defined as the point on the DET 
curve where P mist = P FA , is 4.8. Naturally, more importance should be given to P miss rather than to 
P FA . Therefore, it is assumed that in a practical system the selected activity point on the DET curve 
will be where P m iss=2% and P F a=9%- 
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EXAMPLE 6: DISCUSSION 
From the practical point of view, these examples have illustrated methods and apparatus for 
detection of OLI. An algorithm for detection of OLI by monitoring lungs sounds was developed. In 
order to examine the algorithm performance, a database of recorded breathing sound signals of 
25 patients during OLI and TRI situations was established. It has been shown that assuming a MIMO- 
AR model and selecting the second highest eigenvalue of the residual covariance matrix as a feature 
proves itself as a reliable method for detection of OLI on real breathing sound signals. 

Because of the fact that a pre-processing according to the surgery conditions has to be 
performed, it is disclosed that optional automatic training of the system before every surgery in order 
30 to enable it to set the optimal gain for each microphone is advantageous. 
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Appendix A 



Proof of Equation (9) 

Maximization of (8) with respect to the unknown parameters, AA is achieved via equating the 
sponding partial derivatives to zero. Only the last term of (8) which is: 



corre 



1 N 

n-1 



(14) is relevant to 



l^^-A^W/R-'Cyt^-A^W)]' 

y^]R-M»]-y Mr WA r R-y w) W-/W R " ,A y (w) M+ 

(Ay { ^[n]) r R-'(Ay (A °[»]) 
calculate the derivative of the log-likelihood with respect to A. The derivative of a scalar a with 

respect to a matrix B is defined as: 



3s 
5B 



Bp 

at, 



i 



a*. . . , 



S3 



, where B = 



• « 



b 



UN J 



We shall use the following identities: 



3 [*'A y ]-A[ y 'A'*]"* r 



dA 



A [( Ax) r C(Ax)] = (C + C r X Ax)x r 

If A is a square matrix and/A) is a scalar function, then: 

aiog_rtA)_ 1 d/(A) 

aA /(A) aA 

If A is a square and invertible matrix, then: 

3 A J = |A|(A _, ) r . 



(15) 



(16) 



(17) 



(18) 
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and we obtain: 



dA 



log(/(y[l] y[^]|R.A) = 



(19) 



10 



2LK' 



R 



-i 



J ^ n<0 



y[n)y r [n] - Ay (M) [n]y r [»] - y [n]y (M)T [«] A T 
+Ay (w) [»]y (W)r [«]A r 



(20) 



in order to find the ML estimator of A and R, the above derivatives should be equated to zeto. 
Since R is a covariance matrix then, R = R T and R"' = R r . Therefore we obtain two matrix 
equations with two unknown matrix variables, A and R"' : 

-|i[-R" , yt«ly (w,r t»] - R-*y["]y (A ' )T ["]+ 2R"Ay (w) [»]y (w,r [«]] = o (21) 



2 2tr 



y[% r M - Ay (A ° Wy r M - y[»]y (w)r MA r + 

Ay («) [; , ]y ^)r [/l]A r 



= 0 



(22) 



Eq. (22) can be simplified to: 



nr = ty[n]y r M- A|;y (w) [»]y r [»l-fty[»]y (M)r wl A r 



+ A 



r £y<"W" ,r ["]]A r • 



(23) 



Extraction of A from (21) while assuming the matrix Jy (JW) [n]y (w)ir [«] is invertible, leads to 

n-1 



(9a). Substituting (9a) into the last term of (23), and Extraction of R leads to (9b). 



15 



WO 2005/094179 



PCT/TL2O05/000369 



19 

Appendix B 
Proof of Equation (10) 

Substituting A and R into (8) leads to (24): 

]og[/(y[l] srl*niR.A)]-~ log(2^)- T log|R|- 

-lf[(yW- AuWf R-(y[»]- Aa[n])] ^ ^ ^ 

(24) is a scalar, and therefore the trace operation can be performed on it: 



term in 



J£{gm-a«*»]^ 

\tr=* n 
-1 



■7=1 
ft 



(25) 



The fifth equality in (25) arises by substituting (9b). 

The determinant of matrix is the product of its eigenvalues. Therefore, recalling {/, £ the 
iaenvalues of R in descending order (/,>/,>->/,), |r| can be simplified into: 



w-n* 



1=1 



(26) 



As appears in [16], the smallest L-K eigenvalues of R are used to estimate <x» . Therefore, (26) 



turns into: 

irI ~~ " ^ " 



* ( A V 

-El* lArSH 



(27) 

I. • 

in 



Because of the assumption that 0 > is known, it shall be substituted instead of its estimation 
(27), and therefore (27) turns into: 



WO 2005/094179 



PCT/IL2O05/000369 



20 



*hri'.V) 

(=1 



(28) 



Substituting (25) and (28) into (24) leads to (10). 
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« t l present invention has heen pariicu.ariv desert P-» - * ^ 

■ , ** m anv variations and modifications can be made. Therefore, the ,nvent,o„ ,s not to he 

of -dte invention wil, he m ore tea* u*det*ood hv refetence to the d« fo„ow. 



